Well-localized edge states in two-dimensional topological insulators: ultrathin Bi films 
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We theoretically study the generic behavior of the penetration depth of the edge states in two-dimensional 
quantum spin Hall systems. We found that the momentum- space width of the edge- state dispersion scales with 
the inverse of the penetration depth. As an example of well-localized edge states, we take the Bi(l 1 1) ultrathin 
film. Its edge states are found to extend almost over the whole Brillouin zone. Correspondingly, the bismuth 
(111) 1-bilayer system is proposed to have well-localized edge states in contrast to the HgTe quantum well. 

PACS numbers: 73.43.-f, 72.25.Dc, 73.20At, 85.75.-d 
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Introduction- The quantum spin Hall (QSH) phase 111 121] is 
a new state of matter predicted theoretically, and has received 
a lot of attention recently. This phase is a nonmagnetic insula- 
tor in the bulk or film, and has gapless surface or edge states. 
The edge states consist of counterpropagating states with op- 
posite spins. The notable feature of these edge states is that 
they are topologically protected; they remain gapless even in 
the presence of nonmagnetic impurities and interaction |13|,0]. 
We still know few systems in which the QSH phase is real- 
ized. The first theoretical proposal for the QSH system on the 
Bi ultrathin film by one of the authors |50 . In addition, HgTe 
quantum well has been theoretically proposed jd], and exper- 
imentally shown to be in the 2D QSH state iAm- 

The edge states are localized near the edge, but their pene- 
tration depth £ into the bulk varies between the systems. The 
observation and control of the edge states crucially depends 
on the penetration depth, and it is an important issue how they 
are determined in various systems. In the present paper we 
study the behavior of penetration depth £ in QSH systems. 
From a simple model we show that the minimum penetra- 
tion depth (which is typically reached in the middle of the 
bulk gap) scales with the inverse of the extension of the edge 
states in /c-space. Therefore, if the edge states exist only in a 
small region in /c-space, £ is long. By extending this conclu- 
sion to generic cases, we expect that the penetration depth is 
of the order of the lattice constant, if the edge state extends 
almost over the whole Brillouin zone. To see this, we numer- 
ically study topological properties of bismuth ultrathin films 
and their edge states. Among bismuth thin films, only two 
thin films are proposed to be insulating in the bulk: the (1 1 1) 
single (1) -bilayer film [9] and the {012} 2-monolayer film 
lildl . By using tight-binding Hamiltonians obtained by first- 
principles calculations, we found that (111) 1-bilayer film is 
in the QSH phase and {012} 2-monolayer film is not. We 
also found that the edge states in Bi (111) 1-bilayer film are 
well localized near the edges, compared with the HgTe quan- 
tum well. From these studies we conclude that the penetra- 
tion depth £ corresponds to the inverse of the /c-space width of 
edge- state dispersion. 

Penetration depth of the edge states- We use the Hamilto- 



nian for the HgTe quantum well. 
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where H (k) = ekl2 + d^ (k) a^. Here, I2 is a 2 x 2 unit ma- 
trix, cFa the Pauli matrices, e\^ = C — D {k^. -\- ky),d^ = Akx, 
d^ = Aky, and d^ =M{k) = M-B{kl^ k^). The con- 
stant C is set to zero since it is an overall energy offset. The 
eigenenergies are then given by —Dk'^ ± |d(k) |. Thus D rep- 
resents the asymmetry between the valence and the conduc- 
tion band dispersions . The bulk gap at k = is given by 
2M. In order to consider the edge state on a single edge, we 
deal with a system on a half-plane of ^ < 0. This considerably 
simplifies the results, compared with the ribbon of finite width 
1I11I1 . As edge states only the solutions with e^^ (ReA > 0) 
are allowed. The secular equation 
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where B± = B ± D, gives two allowed values for A: 



A = Ai 



kl^F± y^F^-{M^-E^)/{B^B_), 



(3) 
where F = — ~^)^ ^^ ^ . If we impose a boundary condi- 
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tion tlj{y = 0) = as in llllll , we get 



A1A2 - ^^ 



kt, Ai + A2 = , ^ ^ . (4) 
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From Eqs. ([3]) and Q, we obtain an exact form for the disper- 
sion of edge states 
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The signs correspond to the two branches of edge states with 
opposite spins. Because they are related with each other by 
Kramers theorem, we henceforth consider only the plus sign 
in ([51). Putting ^ into dU), we get 
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, Ai + A2 = 27V, (6) 



where N = A/(2^Bj^B-). These determine Ai^2. If we 
put Ai > A2 , \p ^ gives the physical penetration depth ^ as 
discussed in 11 111] . At the points with A2 = 0, the edge states 
have infinite penetration depth and become bulk states. From 

© this occurs when k:, = k^ = ^{l± ^/iT^^) . 

It can be checked that the states at kx = k^ are located at 
the band edge of the projection of the bulk band, and at these 
points the edge dispersion (|5]) is tangential to the bulk band 
projection. We can rewrite as A1A2 = —{kx — k^){kx — 
k~). Therefore A2(= ^~^) is expressed as ^~^ = A^ — 

y N'^ + {kx — kt){kx — kx). Hence the behavior of l~^ is 
as shown in Fig.[T] It vanishes at the points P± (kx = k^) 
where the edge states are absorbed into the bulk band, and £ 
is minimum when kx = {k^ -\- k~)/2. The minimum value 

Cin is given by ^- J, = N- ^ N^ - {kt - kxf /^- As a 
function of A^, the minimum value of £min is 2/(/c+ — /c~) at 
N = {k^ — k~)/2. This means that the minimum £min of the 
system is roughly given by the inverse of the /c-space exten- 
sion of the edge state dispersion. From Fig. \T\it can be seen 
that the penetration depth £ becomes short when the consid- 
ered edge state is far from the points P± . The inverse of the 
penetration depth £~^ corresponds to an imaginary part of the 
wavenumber perpendicular to the edge direction, and there- 
fore it behaves similarly to the (real) wavenumber. Hence £~^ 
is approximately given by the /c-space distance of wavenum- 
bers from the points P± . 
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HgTe quantum well comes from the fact that the edge states 
are localized within a very narrow region in k space, giving a 
long £. This penetration depth determines the minimal width 
of the system size required for observation of edge states. 

Bi(lll) Ultrathin Film- By extending our theory to generic 
types of edge states, we can expect that the inverse of the 
penetration depth £~^ well scales with the /c-space distance 
from the absorption point P± into the bulk band. Therefore, if 
the edge states extend over the Brillouin zone, the penetration 
depth of the edge states is as short as a few lattice constants. 
We will theoretically show that Bi(l 1 1) ultrathin film is a QSH 
system having edge states with such a short penetration depth. 

For the calculation, we use a tight-binding model con- 
structed from maximally localized Wannier orbitals 1I12I] ob- 
tained from first-principle calculations [13]. The Fermi en- 
ergy lies in the 6p-like states, comprising three conduction 
bands and three valence bands. Therefore, in constructing the 
Wannier orbitals we only retain these six bands. From these 
Wannier orbitals including the lattice relaxation effects of the 
ultrathin films, we construct tight-binding models keeping up 
to third-neighbor hopping amplitudes. 
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FIG. 1 : Penetration depth £ for the effective model with ribbon ge- 
ometry. CB (VB) represents the bulk conduction (valence) band. The 
plot on the right corresponds to a more asymmetric situation, leading 
to a larger £ at the crossing point of the edge states. 



In the HgTe quantum wells, the 2D quantum spin Hall states 
are confirmed by transport measurements (tI |8|] . The pen- 
etration depth of the edge states in these systems has been 
calculated to be relatively long £ ~ 50nm uM- In our 
theory, by plugging the parameters into our results, we get 
£ = 56nm at kx = 0, in agreement with lllin . The coeffi- 
cient D gives rise to an asymmetry between the conduction 
and the valence bands, and the edge state is also asymmet- 
ric: k^ ^ -k-, k^ = 0.62nm-\ k' = -0.024nm-^ 
Thus the penetration depth £ is shortest not Sit kx =0 but 
at kx = {k^ + k-)/2 = O.SOnm-^ with Cin ^ 6.2nm. In 
our interpretation, the relatively long £ of the edge states in 



FIG. 2: (a) Unit cell and latdce vectors, and (b) TRIMs of Bi(lll) 
ultrathin film. The TRIM Fi is the F point, and the three TRIMs 
r2, Fs and F^ are the M points, (c) Bulk energy bands and the par- 
ity at the TRIMs for a Bi(lll) 1-bilayer. (d) and (e): Energy bands 
of the Bi(lll) zigzag and armchair edge ribbons, respectively, with 
a width of 20 unit cells, calculated from tight-binding model, (f) En- 
ergy bands of a eight-unit-cell wide Bi(l 1 1) zigzag edge ribbon from 
first-principles calculations. The size of the symbols corresponds to 
the weight of the states in the edge atoms. 

Figures m (a - c) show the unit cell and lattice vectors, re- 
ciprocal vectors and TRIMs, and the energy band of Bi(lll) 
1-bilayer respectively. Since this system is inversion symmet- 
ric, all states in Fig. [2] (c) are doubly degenerate. This sys- 
tem is proposed to be a nonmagnetic insulator with a bulk 
gap of ~ 0.2eV [9]. We will calculate the Z2 topological 
number u, and if it is nontrivial it is a QSH insulator. The 
wavenumbers satisfying k = — k (mod G) are called time- 



reversal-invariant momenta (TRIM) k = F^ {i = 1,2,3,4). 
For inversion- symmetric systems, the Z2 topological num- 
ber u is defined by (-1)^ = [iti II^^i 6m(r.) (= ±1), 
where ^2m (r^) (= ±1) is the parity eigenvalue of the Kramers 
pairs at Ti and n is the number of the Kramers pairs of eigen- 
states below the Fermi energy lll4ll . The parity eigenvalues 
at the TRIMS Ti{i = 1, 2, 3, 4) are given in Fig. [21 (c) and 
yield the topological number v = 1. We note that both the 
first-principle calculation (without a tight-binding model) and 
the calculation of the Liu- Allen tight-binding model llisll give 
z/ = 1. In Ref. [50 a (1 1 1) 1-bilayer bismuth is proposed to be 
in the QSH phase, from a simple truncation of the 3D tight- 
binding model lllSn . We thus confirmed that the conclusion in 
Ref. ISO remains unaltered in first-principle calculations. 

If we neglect the out-of-plane coordinate, the (111) 1- 
bilayer film has a honeycomb structure. Therefore, as in 
graphene we refer to the two types of simple edge shapes as 
zigzag and armchair edges. Figure [2] (d)(e) shows the energy 
bands of zigzag and armchair edge ribbons of the Bi(l 1 1) 1- 
bilayer. Due to inversion symmetry, all the states are doubly 
degenerate, and they have opposite spins, localized on the op- 
posite edges. In both figures, the number of Kramers pairs of 
edge states on the Fermi energy per one edge is odd, confirm- 
ing that Bi(lll) 1-bilayer is a QSH system. We checked that 
for the zigzag-edge ribbon our result from the tight-binding 
model (Fig. [2] (d)) and that from a first-principle calculation 
(Fig. [21(f)) are in good agreement. 

These edge states are quite different from those in a HgTe 
quantum well, where the edge states exist only near the /c = 
point 161]. Within our calculation, the edge states extend al- 
most all over the whole Brillouin zone. At the Fermi energy 
there are three Kramers pairs of edge states. Thus, the con- 
ductance in a ribbon geometry becomes G = Qe'^/h for a 
clean system. When nonmagnetic disorder is increased, some 
of these edge states become gapped due to elastic scattering, 
while at least one pair of edge states remain gapless, giving 
the conductance ofG = 2e^ /h. These edge states form per- 
fectly conducting channels, similar to those in the graphene 
nanoribbons llldl . In graphene, perfectly conducting channels 
are formed only in the absence of short-ranged disorder; in the 
Bi (1 1 1) 1-bilayer nanoribbon the perfectly conducting chan- 
nel exists irrespective of the nature of nonmagnetic disorder, 
and it gives a universal behavior realizable in experiments. 

Bi{012} ultrathin film- For inversion asymmetric systems 
such as Bi {012} 2-monolayer film, the calculation of v is 
complicated because the phases of the Bloch wavefunctions 
in the entire Brillouin zone are involved 111 livll . The phase 
of the wavefunction is a gauge degree of freedom and can be 
chosen arbitrary for each k. Hence, a simple discretization of 
a formula for continuous k suffers from numerical instability 
due to this gauge choice. 

Hence, we adopt a gauge-invariant discretization method 
proposed in Ref. llisll . It is a merit of the method that we do 
not need to determine the phase of the wavefunction smoothly 
in k space. The mesh size 6kiSk2 should be fine enough to 
satisfy |i^(k)| SkiSk2 < tt at any mesh, where F{k.) is the 



Berry curvature, and 6ki , 6k2 are the width and height of a 
mesh, respectively llisll . This quantity is largest when k is at 
the direct gap k = kg, and the critical size is approximated 
by the k- space nominal size of the band extremum at k = kg. 
From the band structure of Bi{012} 2-monolayer, the criti- 
cal mesh number n^ is estimated to be ~ 100. For various 
mesh numbers exceeding n^ we get the consistent result that 
the Z2 topological number is i^ = 0. Therefore, Bi{012} 2- 
monolayer is an ordinary insulator. 
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FIG. 3: Energy bands of Bi{012}: (a) zigzag and (b) armchair edge 
ribbons with a width of 20 unit cells. 'U' ('L') means that the state is 
localized on the upper (lower) edge. The crystal structures near the 
upper edge are shown in the right panels with 1,2,3 and 4 represent- 
ing the lattice sites. We note that these 4 sites do not lie on the same 
plane. The shaded regions are the bulk energy bands. 

The edge states for ribbons with two types of edges of 
Bi{012} 2-monolayer ribbons and their energy bands are 
shown in Fig. [3l The two edge shapes can be called zigzag 
and armchair edges, although the lattice structure is quite dif- 
ferent from graphene. The number of Kramers pairs of edge 
states on the Fermi energy is even at each edge and this is in 
agreement with our result that the Z2 topological number is 
u = 0. In the armchair-edge ribbon, the edge states are al- 
most degenerate because of the equivalence of two edges of 
the ribbon via mirror- symmetry. These two states have an en- 
ergy difference, due to hybridization of the edge states at the 
opposite edges. Nevertheless, for a ribbon wider than £, the 
energy difference is exponentially small. On the contrary, in 
the zigzag-edge ribbon, the edge states are not degenerate be- 
cause of inequivalence between both edges of the ribbon. 

Penetration depth of the edge states- We now calculate the 
penetration depth I of the edge states of the Bi(l 1 1) 1-bilayer 
film. For the edge states on the zigzag edge of the (111) 1- 
bilayer film, the result is shown in Fig.[4j(i). The penetration 
depths £ of the edge states (including those on the Fermi en- 
ergy) are typically several lattice constants. Hence, for trans- 
port experiments the width of the sample has to be larger than 
a few lattice constants. 

These results on Bi(lll) film agree with our theory on 
I. According to our theory, the penetration depth ^ be- 
comes short when the edge states are distant from the points 



P± where the edge states merge into the bulk (circles in 
Fig-lllii)(iii))- Hence, £ is longer for the states at E^ in Fig.|4] 
(ii), and shorter in Fig.|4](iii). This information is relevant for 
transport which is governed by the states at the Fermi level. In 
Bi(l 1 1), the edge state travels almost over the whole Brillouin 
zone (BZ). Therefore we estimate £ ~ (size of the BZ)~^ ~ 
(lattice spacing), in agreement with the results in Fig.Ui) 

Bi(lll) 1-bilayer film cannot be described by an effective 
model near k = like ([T]). The effective model ([T]) is de- 
rived when the QSH system is described as a band inversion 
between two doubly-degenerate bands, such as HgTe quan- 
tum well, or Bi2Se3. In bismuth ultrathin films, the involved 
bands are Px, Py, Pz orbitals, and the valence and conduc- 
tion bands have different mixing coefficients for these or- 
bitals. Therefore, it is not a mere band inversion, which is 
the reason why the case Fig.lUiii) is realized in bismuth films. 
Bi2Te3 and Bi2Se3 ultrathin films also have edge states sim- 
ilar to Fig. HJiii) [19]. In these cases, however, some edge 
states are close to the bulk bands, leading to very long pene- 
tration depths of about a hundred times the lattice constant . 
We note that our theory assumes isotropy between the direc- 
tion along the edge/surface and that perpendicular to it. For 
layered materials such as Bi2Se3 and Bi2Te3, the penetration 
depth perpendicular to the layer cannnot be predicted from the 
surface-state dispersion in the layer because of the anisotropy. 




FIG. 4: (i) Penetration depth of the edge states on the zigzag edge 
of the Bi(lll). (ii)(iii) Examples of the edge states in the 2D QSH 
systems. 

These short penetration depths of edge states in the Bi (111) 
1-bilayer film are ideal for observation by STM/STS and con- 
trol of the edge states. Furthermore, it is also favorable for 



edge thermoelectric transport [2u\. To utilize the perfectly 
conducting channels of edge states for thermoelectric trans- 
port, short penetration depth is an important factor, because 
longer penetration depth mixes the states at different edges 
for narrow ribbons, and destroys the coherent electron trans- 
port at the edges. 

Conclusion- We analyze a generic behavior of the pene- 
tration depth of the edge states in two-dimensional quantum 
spin Hall systems. We found that momentum- space distance 
between the edge states and the absorption point of the edge 
dispersion into the bulk band roughly gives the inverse of the 
penetration depth. As an example, we calculate the penetra- 
tion depth of the edge states of Bi(l 1 1) 1-bilayer film, which 
we propose to be a QSH insulator. The penetration depth of 
the edge states in Bi(l 1 1) 1-bilayer film is in good agreement 
with our theory. 
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